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We study the molecular mode coupling theory for a liquid of diatomic molecules. The equations 
for the critical tensorial nonergodicity parameters F™,(q) and the critical amplitudes of the j3 - 
relaxation H™/(g) are solved up to a cut off l co = 2 without any further approximations. Here l,m 
are indices of spherical harmonics. Contrary to previous studies, where additional approximations 
were applied, we find in agreement with simulations, that all molecular degrees of freedom vitrify at 
a single temperature T c . The theoretical results for the non ergodicity parameters and the critical 
amplitudes are compared with those from simulations. The qualitative agreement is good for all 
molecular degrees of freedom. To study the influence of the cut off on the non ergodicity parameter, 
we also calculate the non ergodicity parameters for an upper cut off l co = 4. In addition we also 
propose a new method for the calculation of the critical nonergodicity parameter 

PACS numbers: 61.25.Em, 64.70.Pf, 61.43.Fs, 61.20.Ja 



I. INTRODUCTION 



The mode coupling theory (MCT) of the glass transition is by now an important tool to understand experiments 
in and simulations of supercooled liquids ]5J] . For a long time most of the theoretical investigations concentrated on 
simple monoatomic or binary liquids. All universal and even system specific predictions of these investigations could 
be tested on a quantitative level in a system of hard colloids , which is an excellent realization of a hard sphere 
system and in computer simulations for a binary Lennard Jones system . Details of theory and tests for simple 
glass formers can be found in review articles [p|-^2| and articles cited therein. 

Although the theory was originally formulated only for these simple systems, most of the experimental and sim- 
ulation support came from research on much more complex system (e.g. tri-a-naphtylbenzcn e |13| , Orthotcrphcnyl 
(OTP) 0-0, 0ACa(NO 3 ) 2 0.6KNO 3 (CKN) 0-||, Glycerol ||-||, Salol Hj^], toluene |§ and water @@). 
Also most of the experimental methods used, did not measure density correlation functions or their susceptibilities, 
for which the original theory was formulated. Even neutron scattering experiments in systems consisting of molecules 
whose components have different cross sections for neutrons |]l7| , p3| , p8| do not measure the density correlation function 
exclusively, but a mixture of more complicated correlation functions involving molecular degrees of freedom |3l| ] (see 
also for a single linear molecule). Dielectric loss measurements 0,|4j measure directly the correlation function of 
a tensor of rank 1. Depolarized light scattering |l|,|l|,||,|3| , Kerr effect experiments NMR J3|-|38|]and ESR 

p9| (and references therein) measure correlation function of a tensor of rank 2. The mentioned tensorial quantities are 
all related to orientational degrees of freedom (ODOF) , whereas the original theory j4(| only considered translational 
degrees of freedom (TDOF), i.e. the center of mass motion. But of course, when comparing experimental results on 
complex systems with predictions of the MCT for simple liquids it was always reasonable to argue, that there are 
in every experiment couplings to the center of mass motion. E.g. the reorientation of dipoles measured in dielectric 
loss measurements can induce center of mass motion via a translation - rotation coupling. Also the reorientation of 
the polarizability tensor in light scattering measurements is related to a physical rotation of the molecules and will 
therefore be coupled to the center of mass motion of the molecules as well. A slowing down of this motion due to very 
slow structural relaxations can consequently also indirectly be measured in the mentioned experiments. In addition 
it is perfectly justified to perform tests of the universal predictions of MCT in complicated molecular and polymeric 
systems p| for (3 scaling laws and properties of the a relaxation, like time temperature superposition principle and 
wave vector dependent stretching exponents, since the underlying universal features of the bifurcation scenario should 
also remain valid for molecular systems. 

But beyond the universal aspects MCT aims to be a microscopic theory of structural relaxation. This goal is to a 
large extent achieved for simple liquids. There it was possible to obtain quantitative agreement between experiment 
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and theory for the full dynamic range of structural relaxation (i.e. and a - relaxation) 0^,0] . Recently also a theory 
for anomalous high frequency oscillations (Bose peak phenomenon) was formulated within MCT [^2| . The molecular 
mode coupling theory (MMCT) which is under study for a few years now, intends to extend this line of research 
to experimentally relevant molecular systems. There are three different mode coupling theories for the description 
of different aspects of molecular degrees of freedom. In J3^,[l3| the motion of a single linear molecule in a liquid of 
spherical atoms is studied. In [Q a site - site description is formulated, which is perfectly adapted to study neutron 
scattering experiments of molecular systems. In this approach the atomic structure of the molecules is considered. 
Finally, the MMCT (45), we are using in this work [Q, is devised to investigate the dynamics of a liquid of linear 
molecules. For this purpose a self consistent mode coupling theory for the dynamic correlation functions of tensorial 
densities pi m (q,t) was developed. These densities are the generalized Fourier components of the microscopic density 
p(x, fi) in an expansion in spherical harmonics with respect to the orientation f2 = [9, <p) of the molecules and plane 
waves with respect to x. An extension to arbitrary molecules is given in ]4^| . A theory for arbitrary molecules was 
also formulated in [Q. First results for the tensorial nonergodicity parameters (NEP) and critical amplitudes were 
obtained for dipolar hard spheres Q and [|^| , respectively. A study of the phase diagram for glass transitions of a 
liquid of hard ellipsoids was performed in ref. |30| . Several aspects of the theory for general molecules were tested 
against simulation for water in ref. [5l|j52] | and treating water as a linear molecule in ref. 

As in the MCT for simple liquids the static structure factors S™,(q) in the q frame, i.e. the coordinate system in 
which the z axis is along the wavevector, completely determines the longtime dynamics and thus the NEP and critical 
amplitudes. Note that the structure factors (static and dynamic) are diagonal in m in the q -frame The static 
structure factors have to be known to solve the equations of MMCT. They are either obtained by analytical theories 
as e.g. in ref. [|53|-f55| for ellipsoids or they have to be taken from simulations. In this work we present a detailed 
comparison of MMCT calculations of the NEP and critical amplitudes with the results of simulation for a liquid of 
diatomic molecules. A description of the system and the simulation can be found in |5^,|5^]. There also tests of the 
universal properties of MCT are presented. For our comparison between simulation and theory the static structure 
factors are taken from the simulation. In a preliminary study only the diagonal static correlators S" L (q) were used as 
input and also the dynamical correlators S™(q,t) and thus the NEP were assumed to be diagonal |pq| . This severe 
approximation has lead to unphysical results like the existence of two different transition temperature for ODOF and 
TDOF. In our study we use diagonal and non diagonal elements of the static structure factor as input to calculate all 
components of the NEP and, in addition, of the critical amplitudes. For the calculation of the NEP, we also extend 
the necessary upper cut off l co for the index I to l co = 4. 

The paper is orga nized as follows. In chapter (||) we revi ew the main equations and concepts for the calculation 
of the critical NEP ( II A ) and the critical amplitudes ( II B ). In chapter (111) we discuss the influence of different 
approximation schemes on the theoretically obtained critical te mper atures (III A) and its relati on to simulation results. 
Then we present the comparison of theoretical critical NEP ( HIB ) and critical amplitudes ( [II C ) with simulations. 
Conclusions are presented in chapter IV and an appendix describes how the critical NEP can be obtained from the 
liquid side close the the ideal glass transition. 



II. MOLECULAR MODE COUPLING THEORY 



A. Non ergodicity parameter 

The derivation of the equation of MMCT for the dynamics of linear molecules and general molecules can be found 
in |4(| and |l7j respectively. We only repeat the basic definitions and equations and refer the reader for details to the 
literature. For the present work, in which we want to calculate the critical NEP and the critical amplitudes of diatomic 
molecules, we only need the equations for linear molecules in the limit of time to infinity. The basic quantities are 
the correlation functions of tensorial densities pi m (q) and tensorial current densities jf^fat) 

N 

p lm (q,t) = V^i l J2 e i? ' Snit) Ylm(ti n (t)) ■ (1) 

N 

J,« (£ t) = yfai 1 J2 j* Sn{t) Ylm(tin(t)) ■ (2) 

71 = 1 

The Yi rn {VL n {t)) are the standard spherical harmonics and we follow in our notation the text book by Gray and 
Gubbins |5^] . v% is either the center of mass velocity v n (t) of the n th molecule or its angular velocity uo n (t) depending 
on the index a: 
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where T and R stands for translational and rotational part, respectively. For the calculation of the NEP we in 
principal need all spatial components of the currents but here we use as in p6| only the projection on directions 
defined by the wave vector q and the angular momentum operator L. Taking into account also transversal currents 
will lead only to small correction in the NEP (60|. We therefore define the longitudinal currents jf m (q,t). 

jLm) = ^(q a J a ) i (<f) , ae{T,R} (4) 

Qi \ / Irn 



with 



and the definition 



q , a = T 

L , a = R (6) 

. The quantities we are going to calculate are NEP FJJt(q) in the q -frame i.e. in a coordinate system in which the 
z-axis is given by the direction of the wave vector q = (0,0,(7). in this coordinate system all correlation functions 
{Pi m (q,t)pi' m '(q,0)) — 5 rnm 'S" 1 l ,(q,t) are diagonal in m. They are real and depend on \m\ only. The same holds for 
all other tensorial quantities, we will use in the q - frame. The NEP are given by 

F$(q)= Km S%(q,t). (7) 

t — >oo 

As input for the mode coupling equations we need the static structure factors S™,(q) They are directly taken from 
the simulation of Kammerer et. al. pq . The off diagonal elements of the static structure factors were not yet published 
and had to be determined from the raw data. 

The equation for the matrix F(g) = (F{p(q)) of NEP can then be written as 

F(q) = [l + K(q)S- 1 (q)]- 1 S(q) . (8) 

The matrix K, is related to the Laplace transform of the dynamic current correlation function and can be expressed 
as the inverse of a memory matrix J^^q) = (T™\qj) at t = 00, 



t\m = xyfa) (^ m (<r% Q qP' (<?)■ (9) 

OLOi' 

The mode coupling approximations yield: 

""'<«> -sKs)' 

x Vwtu ll2 v 2 (Q, 92; m, m 1: m 2 ) F™} F™? (q 2 ) , (10) 

Vm^i^iq, qi, 92; m, m 1 ,m 2 ) := 

Vu x i 2 {q,qi,q2;m,mi,m2)-vp i ,j, 2 (q,q 1 ,q 2 ;m,mi,rri2)* , (11) 

Vu 1 i 2 (q,qi,q2;m,mi,m 2 ) := 

E <>h («« 9i,Q2; m, mi ,m 2 ) cH ( gi ) + (-l) m (l <— 2) (12) 

h 



with 



and 
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u?i 1 i 2 (q,qi,q2;m,m 1 ,m 2 ) 



(2/! + l)(2/ 2 + l) 
(21 + 1) 



i[l + (-l)^] 
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Here the functions C(l\l2l;mx,m2,m) are the Clebsch Gordan coefficients and d l m , m (Q) are related to Wigner's 
rotation matrices (we follow the notation of Gray and Gubbins). For given Euler angles = ($, O, x) they are defined 
as H) 
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9i is the angle between <f and <fi. The prime at the first summation in Eq.(|10|) restricts q\, q\ such, that q\-\-q\ = 9- 
Eqs. (|| - and ([h]) form a closed set of infinitely many coupled nonlinear equations for the NEP. To obtain a solvable 
theory we have to restrict I to be smaller than an upper cut off, I < l co . The resulting equations can in principle be 
solved by a fixed point iteration algorithm. The physical control parameter like the temperature and the density only 
enter via the static structure factor. At a critical temperature or density the solution of this equations bifurcates from 
all functions F^{q) being zero to nonzero. In the simulations of Kammerer et al. |56|,|57j the temperature was used as 
control parameter. Close to the transition temperature T c the stability matrix of the iteration (see below) will have a 
largest eigenvalue Eq approaching Eq = 1 from below. Consequently the convergence of the iteration is getting very 
slow close to T c . The time for one iteration depends very sensitive on the upper cut off l co . For l co — 2 one iteration 
took 10 minutes on a MIPS R10000, for l co — 4 this time increased to 6 hours. We therefore concentrated on l co = 2, 
to determine the transition point with a high accuracy and used the calculation for l co = 4 mostly as check for the 
sensitivity of our results against changing the cut off. 

To overcome some of the restrictions connected to the critical slowing down of the convergence close to T c we 
determined the critical NEP, i.e. the NEP at T c , with two alternative methods. For the standard fixed point iteration 
we started at a temperature low enough to be in the glass state. Then the temperature is increased very slowly. At 
every temperature the equations for the NEP are solved by the iteration 



F^ n+1 \q) = g(F n ,e) 



(15) 



where g(F( n >) is the right hand side of Eq. (|) and e = (T c - T)/T c . This iteration converges exponentially fast 
towards its solution as long as the temperature is not the critical temperature. The convergence rate is determined by 
the largest eigenvalue Eq of the stability tensor C\y = dQ\/dF\i [plf . The index A is an abbreviation for wavevectors 
q and rotational indices /, I', m. The exponential convergence rate is then lni^o- Close to and below T c the eigenvalue 
Eq can be written as Eq = 1 — Ayfi with A being a positive constant. Therefore the convergence rate is Ayfi and 
the number of iterations to obtain convergence diverges inversely proportional to y/e close to T c . With this number 
of iterations the deviations of our NEP from the true critical NEP are proportional to y/e, since the NEP exhibit the 
well known square root singularity (cf. ^jj). 

If the temperature is increased above T c , i.e. e < 0, there is no non zero solution for the iteration Eq. (|15|). 
Nevertheless for < — e <C 1, the iteration is nearly stationary for a large number of iteration of the order |e| / 2 
(see appendix). The approximate critical NEP is determined as the stationary point F(e) whose change along the 
eigenvector with eigenvalue 1 of the critical stability matrix is minimal during iteration. But contrary to the NEP 
determined from the fixed point iteration for T < T c , the stationary NEP differ only in order |e| instead in order 
yfi from the true critical NEP. Consequently this property allows us in the following to crosscheck the very accurate 
results for l co = 2 obtained from the fixed point iteration and also to obtain the critical NEP for l co = 4. 



B. Critical amplitudes 



A central prediction of the mode coupling theory of the glass transition in simple liquids is the existence of the 
(3 - relaxation regime |^,^,^|. For the problem of a single dumbbell in an isotropic hard sphere system [ j32| it was 
demonstrated, that the /3 relaxation law can be detected in every quantity, which couples to the density. For a liquid of 
anisotropic molecules it is not very well defined which degree of freedom is driving the glass transition. The equations 
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of MMCT couple all degrees of freedom and there are situations where the transition is not caused by the TDOF, but 
the ODOF |3CJ. But even in these systems the factorization theorem is generically valid for all correlators. This can 
be proven using the standard techniques ||. Therefore every dynamic structure factor S(ft(q, t; T) can for — 1 <C e <C 1 
in the f3 relaxation regime be written as 

S%(q,t;T) = FlP(q) + H% (q)G(t/t ; a) (16) 

The function G(t/to) is the same for all l,l',m,q. to is an overall microscopic scale. H\J L ,{q) are the critical 
amplitudes determining the intensity of the asymptotic f3 - relaxation for a certain combination of I, V , m and q. Also 
the correction to the asymptotics, which determine, besides the temperature, the range of validity of the law Eq. jl^ ) 
depend on these amplitudes (cf. ref |34| for simple liquids). Differences in the observability of the critical correlators 
between depolarized light scattering experiments and dielectric loss measurements, can be explained by differences in 



the amplitudes H[[i(q) involving I = 2 and I — 1 respectively (see ref. |43 for a single molecule). To determine the 
amplitudes numerically to a high precision it is necessary to be very close to the transition point, to make sure that 
all correction terms of order e are small compared to the leading term of order i/e. Due to the difficulties described 
above we could only determine the critical amplitudes for upper cut off l co = 2. 



III. RESULTS 



A. The critical temperatures 

For l co = 2 and in the full diagonalization approximation H|, in which the static structure factors S(q), the glass 
form factors F(q) and the memory matrix T(q) are assumed to be diagonal with respect to I, the glass transition 
temperature for the TDOF predicted by MMCT is below the transition temperature of the MD simulations T C M£> = 
0.477. Note, that this temperatures are given in Lennard Jones units (cf. |36|]). In all other known examples the MCT 
overestimates the tendency for vitrification. As an additional artefact of the full diagonalization the ODOF vitrify at 
a lower temperature than the TDOF. Since the top down symmetry of the dumbbells is broken, the full equations of 
MMCT (^) - ([To]) do generically not allow for such a scenario. As soon as we take S(q) non diagonal all degrees of 
freedom undergo a glass transition at the same temperature above the MD result. To study the influence of different 
diagonalization approximations a bit more in detail, we investigated several cases, with the main condition of S(q) 
being non diagonal: 

1. F(q) and F(q) diagonal (dd) 

2. F(q) diagonal and F(q) non diagonal (dnd) 

3. F(q) and F{q) non diagonal (ndnd) 

Let us discuss l co = 2 first. For this case T c has been determined very accurately from the asymptotic behavior 
(1 — Eq(T)) 2 cx e (see chapter II B) for the largest eigenvalue Eq(T). Fig. |l|. demonstrates this law for case <^). 



The highest transition temperature is obtained for case (|lj). Here the transition temperature is roughly three times 
as large as the MD result, T dd = 1.4. In case (||) it only slightly decreases to T dnd = 1.38. If everything is taken 
non diagonal (case [|) the transition temperature is T c = 0.7521. Although still twice as large as the MD result, the 
discrepancy of this result from T^ ID is comparable to other known cases and consistent with the usual 20% accuracy 
of the critical density jj| . The equations are too complex to get a deeper theoretical understanding of this seemingly 
erratic jumping of the transition temperature, dependent on the approximation we are using. Particularly the fact, 
that the vertices (Eq. |i"2| ) are not positive anymore, makes an analytical prediction impossible. But it is at least 
possible to rationalize the behavior using a combination of physical and mathematical arguments. First of all it is 
quite clear that the full diagonalization where all matrices are assumed to be diagonal is too crude to describe the 
coupling of TDOF and ODOF for the system of diatomic molecules. The TDOF and ODOF are only coupled via the 
diagonal memory function aa (q). The coupling of the equations for different I is considerably reduced compared 
to the case where S(q) is taken to be non diagonal. E.g. in the ^oo TT (9) component of the memory matrix only 
terms of the form (symbolically) J2u m ' V^vi> ) 2 appear in the full diagonalization approximation, since the Clebsch 
Gordan coefficients C(0l'l", 0, to', to"), which enter into the vertices (cf. Eqs. ( p^| ) - (fl3|)) are nonzero for I' — I" only. 
Similarly, the memory functional J-"l aa (q) only contains couplings of the form F^F]^ and F^F^. For J 7 ^ aOL (q) 
the Clebsch Gordan coefficients in the vertex allow "self" couplings (F^) 2 , (F^) 2 and couplings to I = in the form 
FqqFqz, but no "self" couplings (-Foo) 2 - Due to the absence of (Fg ) 2 in T"^ aa (q) a freezing of the center of mass 
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motion, i.e. I = V = does not imply a freezing for quadrupolar dynamics I = I' = 2. If the vertex for the coupling of 
the NEP with / — 2 and Z = in and are not large enough, exactly this structure of the memory matrix allows 
generically a separate transition of the I = and the I ^ components of the diagonalized dynamic structure factor as 
observed in f58fl . But we have to stress that the approximation is not per se inadequate. In the case of water M] the 
full diagonalization approximation leads to a rather satisfactory agreement with simulations, without the artefact of 
separate transitions and too low transition temperatures. The pronounced angular dependence of interaction between 
water molecules, which is reflected in the fact that the static structure factors for I — 2,m = 0, 1,2 and I — are 
of the same order, yield large enough vertices to produce a single transition temperature of the TDOF and ODOF. 
In the present case, the structure factor Sq is clearly more dominant than (see Fig. |J). This is different from 
water where all m for I = I' = 2 are important. This leads to the a posteriori conclusion that in general the full 
diagonalization can only be used (if at all) for systems with "very strong" static translation rotation coupling. This 
statement, unfortunately, cannot be further quantified. 

If we now take the non diagonality of S(q) serious but leave all other matrices diagonal (case |l|), additional coupling 
between TDOF and ODOF appear, which may lead to an effectively stronger coupling. Although the equations for 
the different / components of the NEP still couple only via the diagonal memory functions FJJ 1 °" (<?); the diagonalized 
memory matrix contains now static couplings between all NEP. E.g., aa (q) contains additional couplings between 
the TDOF - correlator Fq and the correlators involving I = 2 and, even more important, " self coupling" terms (Fq ) 2 
due to the non vanishing structure factors Si and S^(q), respectively. This of course does not explain, but at least 
makes plausible the dramatic increase of the transition temperature. Any slowing down of let's say the TDOF is 
immediately transferred to all other degrees of freedom and causes a further slowing down of the TDOF due to the 
feed back via the memory function. This enhances the tendency towards vitrification and also is responsible for the 
existence of a single transition temperature. 

The reason for the decrease of the transition temperature, when we give up the diagonalization approximations 
for J-{q) and F(g) (case ||), is not obvious. We only note, that the transition temperature is decreasing from case [j] 
(T c = 1.4) over g (T c = 1.38) to || (T c = 0.752) i.e. the more off diagonal elements of the matrix T are taken into 
account. J- is the t — > oo limit of the memory matrix i.e. the random force correlation function. Therforc it seems, 
that the more components of the random forces are coupled, the lower is the transition temperature. This implies, 
that the more the random forces can mutually influence each other, the more difficult it is to form a glass. Although 
we cannot proof this statement on mathematical grounds, it describes an feasible physical phenomenon. 

To test the sensitivity of our results to changes in the cut off, we also solved the MMCT for upper cut off l co = 4. 
The larger cut off value for I reduces the transition temperature further towards the simulation result. For l co = 4, 
an upper and a lower bound for T c has been determined. The lower bound is the highest temperature for which the 
NEP are still nonzero after 88 iterations. The upper bound is the temperature for which the NEP are converging to 
zero after about 24 iterations. Since the time per iteration increases dramatically upon increasing the upper cut off 
the transition temperature could only be determined within 5%, T c = 0.61. It is encouraging that the real transition 
is approached upon increasing the cut off, but our arguments presented above, show, that this is not necessarily the 
case. Which of the competing mechanisms influencing the transition temperature is dominant, cannot be predicted 
on general grounds. 

B. The non ergodicity parameters 



In the following we concentrate on the results for the normalized NEP f"J(q) = F^ {<£)/ \J Sff (q) S^, (q) without 
any diagonalization approximation. Fig. || shows the normalized diagonal terms of the matrix of NEP /™(g) for 
(Z,m) = (0,0), (1,0), (2,0), (2, 1). Not shown are the results for (I, m) = (1, 1), (2,2), since they do not exhibit very 
much structure. The corresponding simulation result is taken from |57fl . It was obtained by fitting a von Schwcidler 
law plus corrections Sff,(q,t) = Fjp(q) - H™{q){t/T a ) b + (H^)^,(q){t/T a ) 2b to the simulation results for the time 
dependent density correlation function S™,(q,t), where r Q is the a - relaxation scale. There are three different 
theoretical curves. The two curves at temperatures T c — 0.7522,0.7521 are obtained with the fixed point method (on 
the glass side of the transition) and the quasi stability criterion (on the liquid side), respectively, as described above 
for upper cut off l co — 2. Their good agreement demonstrates the high accuracy of the solution. The third theoretical 
curve shows the result for upper cut off l co — 4 using the more accurate quasi stability criterion. Compared to the 
results in p8[ a clear improvement of the agreement with simulations can be observed. Especially the q - dependence 
of the functions are very well reproduced. Even a feature like the prepeak in fg (q) at q ~ 2.5 is reproduced as a 
shoulder in the corresponding theoretical result. This peak is not present in the static structure factors. Since Sii(q) 
has a peak at about q ~ 2.5 it could appear due to a dynamic coupling of the ODOF, especially the one involving 
I = 1 and the TDOF. Note also, that the mentioned peak exactly corresponds to the first peak in fii(q). There is a 
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tendency, that the agreement is best around the wave vector, where the structure factor Soo(q) has its first peak and 
is getting worse for large wave vectors. This might be interpreted as an indication for the glass transition being driven 
also for the investigated system of diatomic molecules by the TDOF. From investigations of other systems 
we know that different scenarios are possible. 

Similar to increasing q, the agreement between simulation and theory gets worse with increasing I. This is expected 
due to two different reasons. First higher I correspond to a higher angular resolution and are therefore probably much 
more affected by the mode coupling approximation. Second, higher I are of course much more sensitive to the cut 
off l co than lower I. The curves for larger cut off increase the quality of the comparison with the MD results. But it 
is important to note that in our case no general rule can be given of how much the quality of the results for lower 
values of I can be improved by increasing the upper cut off, as this was done in [3j| . In the case of a single dumbbell 
in a liquid of hard spheres the glass transition temperature is completely determined by the hard sphere liquid and 
does not change by increasing the cut off l co . As explained above, in our case, T c can depend very sensitively on l co . 
But this influences directly the amplitude of the NEP via the trivial effect of the temperature on the static structure 
factors. We already compensate as much as possible for this mechanism by presenting only the normalized NEP. But 
as in the case of hard spheres, there is still the effect, that also the normalized NEP are proportional to the static 
structure factor. This is a very nontrivial phenomenon, since the existence of negative vertices in the mode coupling 
functional J-{q) could in principle lead to a violation of this correlation. But as can be inferred from Fig. || the NEP 
/oo (?) i $21 (?) f° r ^co = 4 are systematically larger than the one for l co — 2 in a large region around the first peak of 
the structure factor 5q ((?) , without big differences in the functional form. This effect, especially for F® (q) where the 
mentioned trend is valid for all wavevectors, can be mainly understood as a consequence of the transition temperature 
being smaller for l co = 4 than for l co = 2, which causes the the first peak of SQ (q) to increase. Additional evidence 
for this reasoning is presented in Fig. [|. In this figure we show the results for l co = 4 at temperature T = 0.60, 
obtained with fixed point method, and T = 0.63, obtained with the quasi stability criterion. The lower temperature 
is still in the glass (The required accuracy i.e. the maximum difference between consecutive iterates is smaller than 
10~ 6 , is reached after 88 iterations). At the higher temperature the accuracy is first increasing as expected (see 
appendix). But after 24 iterations it begins slowly to decrease and after 42 iterations the iterates start to converge 
quickly towards the solution F(q) = 0. The results for the lower temperature agree, except for I = 0, much better 
with the simulation, than the one for the higher temperature. But since the deviations to the true critical NEP at 
roughly the same number of iterations are of order \T — T c \ for the higher temperature compared to order \\/\T — T c \ 
for the lower one, we have to conclude, that the results for the higher temperature are closer to the critical NEP of 
the theory. The better agreement with simulations of the NEP at T = 0.600 is a trivial consequence of the fact, that 
positive NEP increase with decreasing temperature. Due to this influence of the value of T c on the amplitude of the 
normalized critical NEP, we cannot in general conclude that increasing the cut off l co leads to a better agreement with 
the simulation. It might even happen, that the agreement with simulations gets worse instead of better, if increasing 
the cut off would lead to a larger transition temperature. This is possible due to the existence of negative vertices in 
the mode coupling functional T. 

The observed trends allow the reasonable hypothesis, that the temperature effect could be the main source for the 
deviations between simulation and theory. In general the main structural features in the normalized non ergodicity 
parameter are very well represented, but they are systematically to small for nearly all wave vectors, exactly as 
expected, if the theoretical transition temperature is too large. 

Fig. H demonstrates, that the theory also gives good results for the off diagonal NEP. We found that Sq 2 and Fq 2 
are the only important off diagonal components of the static structure factor matrix S(q) and NEP matrix F(q), 
respectively. In Fig. || we therefore show the normalized NEP /q 2 . The quality of the result is even better, than for 
the diagonal components of the NEP. 



C. Critical amplitudes 

The critical amplitudes are determined only up to an overall scale factor. I.e. our theoretical results cannot be 
directly compared to the simulation results. But once we have chosen a scale factor for e.g. the amplitude Hq , all 
other amplitudes should be multiplied with the same scale factor to compare with the simulations. In fig 0(a) we have 
chosen a scale factor of 200 to obtain best agreement with the normalized critical amplitude /in n (<z) = Hnoil) / Son(q) ■ 
The features of this component are the same as in simple glass forming systems [ 
position of the first peak of SQ (q). Simulations and theory compare quite well for 2 < q < 7 and show deviations at 
other wave vectors. With the chosen scale factors the other diagonal elements of the critical amplitude matrix show 
strong deviations from the simulation results. Especially h^q) (see Fig.^(b)) does not even disagree in amplitude but 
also in the form, except for the minimum at q ~ 3. This is not unexpected, since the simulations have shown strong 
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differences between the form of the dynamic correlators involving odd and even I |57|]. Due to the weak top down 
anisotropy of our diatomic molecules, the dynamic of the correlators involving odd I is only weakly coupled to the 

|. 180° jumps are still possible on a much faster time scale than the translational 
motion |37| . Consequently there are strong corrections to the asymptotic results for the correlator (q, t) and the 
amplitude h ( l 1 (q) is not very well defined. 

The deviations between simulation and theory in (see Fig. §(c)) are more serious, since S^iq, t) exhibits a well 
defined (3 - relaxation regime. We can improve the agreement between simulation and theory by choosing a free scale 
factor for the simulation curves. The result is shown in Fig. ^|(d) to demonstrate, in contrast to h\^, that essential 
structural features in h° 2 are indeed reproduced by the theory. As argued above the dynamic correlators and therefore 
also the critical amplitudes involving I = 2 are much more affected by the cut off l co = 2 as the one with lower I. 
This might be the reason for the rather large discrepancy found for h%^. To determine the critical amplitudes, it 
is necessary to be very close to the critical point. Restriction in computer time didn't allow us to determine T c for 
l co = 4 with high enough accuracy to get reliable results for the critical amplitudes. 

The critical amplitudes with m > do not exhibit very much structure. In Fig. |^ we show for completeness the 
result for h\ 2 {q), which could in principle be measured in light scattering experiments. Again we choose an overall 
amplitude prefactor as a free fit parameter, but the agreement is still not very good. Much better is the agreement 
(similar to the NEP) for /iq 2 (see Fig. ||), although we still had to choose the amplitude scale of the simulation as a 
free fit parameter. 

IV. CONCLUSIONS 

We have performed a quantitative test of MMCT for a liquid of diatomic molecules. The static structure factors 
from simulations were used as input for the MMCT to calculate the critical temperature T c , the matrix of critical 
NEP F[p(q) and the matrix of critical (3 - relaxation amplitudes HJfi(q). Several approximation schemes were used 
to test the sensitivity of the results against changing the degree of diagonalization of the F^}{q) in / and V and the 
dependence on the upper cut off l co . As maximum cut off, l co = 4 was used. Since the computational effort increases 
strongly with l COl we used a new, more accurate, method to determine the critical NEP from the liquid side of the 
transition. 

As expected, by giving up any diagonalization approximation for 5 ; ™(q), F{p(q) and the memory functional 
TVu aa (q) all degrees of freedom vitrify at a single temperature. In fact, to obtain a unique transition tempera- 
ture, it is sufficient to keep only the static structure factor SJp(q) nondiagonal in I and I'. The strongest effect of 
successively applying diagonalization approximations for the NEP FJJ} (q) and the memory functional T]$ aa (q) is the 
change of the transition temperatures T c . Also using different cut offs l co changes the transition temperature. 

Contrary to T c , the overall form of NEP is much less sensitive to the different approximation schemes. The 
agreement with simulation is qualitatively quite good for cut off Z co =2 as well as for Z co =4. In some cases especially 
for the I = 0, I' = component and the off diagonal I = 0, /' = 2 component of the NEP the agreement is even 
quantitative for wavevectors around the first peak of the structure factor 5 , Q (q). The comparison of the NEP with 
simulations for the correlator with (/ = 0, I' = 0) is clearly improved using the cut off l co = 4 instead of l c0 = 2. Other 
correlators are better represented for l co = 2. We therefore cannot conclude that in general a further increase of the 
cut off will lead necessarily to still better agreement with simulations. 

Also the wave vector dependence of the normalized critical /3-relaxation amplitudes agrees quite well with the 
one obtained from simulations. Contrary to the prediction of MMCT no common amplitude scale for the critical 
amplitudes with different I and V could be found. At present it is not clear whether this failure is real or due to the 
fitting procedure used in |37| to obtain the beta relaxation amplitudes. Also the restriction to cut off l co = 2 for the 
determination of the critical amplitudes could cause the found discrepancies, since also correlators with / > 2 will 
contribute to the critical amplitudes at I = 2. By neglecting them an error in the amplitude scale is possible. 

In summarizing our study, we may say, that the MMCT offers an overall consistent description of the glass transition 
in molecular liquid of diatomic molecules, at least concerning the critical NEP and the critical amplitudes. 

APPENDIX A: DETERMINATION OF THE NEP USING QUASI STABILITY 

If the temperature is chosen above the critical temperature, the only stable fixed point of the iteration ( |l5| ) is F = 0. 
But there is still the possibility of determining the critical NEP with even higher accuracy than with the converging 
iteration for T < T c . To implement the method, we have to initialize the iteration with a Fjf \ which is close to F^, 
where fi is a superindex for q, I, V and m. This can be achieved by first using the fixed point iteration 
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^({^},e 



(Al) 



below T c and increasing the temperatures in small steps until the liquid regime is reached. The critical temperature 
will eventually be missed by a small value < -e « 1 and the iteration converges after a long transient quickly to zero. 
Since the iteration at a given temperature is always initialized with the final result of the previous temperature, the 
above assumption is fulfilled. The long transient is caused by the marginal stability of the critical point. The iteration 
will be dominated for a long time by the critical direction of the iteration, which is given by the eigenvector e c , which 
belongs to the eigenvalue Eq = 1 of the critical stability matrix (C° /) = {dQ^/dF^Y 61 1. The NEP i^(e) at which 



the iteration is almost stationary is given by minimizing the distance of the projection of AC^ = G ^({F^i , e}) — 



on the critical direction e c 
be written as 



where for technical reasons the left critical eigenvector e c was chosen. This condition can 



dFu.> 



^e^({F M 4,e)-e^=0 



(A2) 



Eq. (A2) determines the quasi stationary point F^e) such, that matrix (C MAl / {{F^n , e)) has an eigenvalue Eo(e) = 1 
for e < 0. e < indicates that the temperature is T = T c (l — e) > T c . This stationary point F^ differs from F^ in 
order e! This can be proven by using an argument analogous to deriving the (3 - relaxation equation [[|. Defining 8^ 
by F^ = F^ + dfj, and using the property that C M/ y ({F^r}, e) is a smooth differentiable function of e, Eq. ( A2) can be 
written, in leading order in 5^ and e as 



° = E^(E c ^v'(f c ,o)v + 



de 



(A3) 



here we have in addition used that is the left eigenvector of the stability matrix at the critical point with eigenvalue 
Eq = 1. The tensor C^i^n is the partial derivative of the stability matrix with respect to F^i . Eq. (A3) immediately 
shows, that <5 M is (at least) of order e). Note, that the argument crucially depends on the well known phenomenon, 
that the static quantities like structure factors, which determine completely the stability matrix {C^i ({F^i>}, e) are 
varying smoothly across the glass transition. 

We now want to estimate the number of iterations required to approach the quasi stationary point F. For that 

(n) 

matter the iterates are expanded around the quasi stationary point 



pin) _ p 

r fj. — A* 



(5, ( , n) 



Neglecting terms of order (S^) 3 the iteration (Al) can be rewritten as 



6 (n+l) = A ^ 



f(™)ir(") 



(A4) 



where the quantities with" are taken at ({F^n},e). The slowing down of the iteration is due to the component o f 6^ 
along the critical direction e c . This component can be extracted by writing 8^ n > — cS n ^e c and multiplying Eq. (A4) 
with the left eigenvector e c of the stability matrix C. The resulting equation for is then 



a 



(n+l) 



+ (A-l)(oW) 2 



(A5) 



Here we made use of (A.2) and have normalized e c ,e c such that J"^, e 
A I V 



a 



1. The quantities a 



£„e£AC„ and 



n" e u,' e t. 



can be identified as the separation parameter and the exponent parameter < A < 1 of 



the (3 - relaxation theory, respectively PJ6q] ■ Since consecutive iterates are very close to each other in the vicinity of 



F we can rewrite (A5) as a differential equation. 



- a = (|<t| + (1 - A)a 2 ) 
Lets assume we start the iteration at t = with a(0) = oq. Then the solution of (|A6|) is 



(A6) 



a(t) 



tan(v'(l - A)|ct| t 



arctan(ao — )) 



(A7) 
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Since the iteration is initialized with a value for F slightly in the glass, the initial deviation ao from F is of order 
y/j^j = 0(y/\e\). From Eq. ( 
F (i.e. a(t) = 0), and will stay close to F (i.e. a(t) < O(y/\o\) for the same amount of iterations before the iterates 
decay to zero. This proves our statement that the quasi stationary solution of the iteration for T > T c is approached 
with the same amount of iterations as the stable solution slightly below T c . But there the glass solution agrees with 
F c only up to order i/e the quasistationary solution F however agrees with F c up to order |e| -C yjef. 



A7) then follows, that we need a number of iteration of the order l/-y/|e| to approach 
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Here we always restrict ourselves to the idealized mode coupling theory, which neglects all ergodicity restoring processes 
like e.g. hopping. 
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FIG. 1. The square deviation of the largest Eigenvalue of the stability matrix (see text) from its critical value Eq(T c ) — 1 



as function of temperature. The full line is a linear fit (1 
magnification for the same quantity very close to T c . 
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FIG. 2. A selection of structure factors S™,(q) of a liquid of diatomic molecules in the q -frame at temperature T = 0.7521 
and pressure p — 1 in Lennard Jones units. The full curve is the structure factor So (q), the long dashed curve is £§2(5), the 
dotted curve is £22(9), tne dashed dotted curve is Si 2 (g) and the short dashed curve is So2(<z) 
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FIG. 3. Normalized diagonal elements fu'(q) of the non ergodicity parameter matrix. As normalization the corresponding 
diagonal elements of the static structure factor are used. Shown are the elements for (l,m) = (0,0) (a) ,(l,m) = (1,0) (b), 
(l,m) = (2,0) (c) and (l,m) = (2, 1) (d). The full line is the result of the simulation at T = 0.477 obtained by fitting a von 
Schweidler law plus corrections of order (t/r a ) 2b to the simulation data. The estimated critical temperature is T c — 0.475. The 
long dashed line is the result of the fixed point iteration for upper cut off l co — 2 at T = 0.7521 slightly below the theoretical 
T c — 0.75215. The short dashed curve is the result of the quasi stationary criterion for upper cut off l co = 2 at T = 0.7522 
slightly above the theoretical T c . The dashed dotted curve is the result of the quasi stationary criterion for upper cut off l co = 4 
at T = 0.630 (above the theoretical T c ) for l co = 4. 
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FIG. 4. Normalized diagonal elements fjj{q) of the non ergodicity parameter matrix for l co = 4 obtained with fixed point 
method and stationary criterion respectively in comparison with simulation. As normalization the corresponding diagonal 
elements of the static structure factor are used. Shown are the elements for (I, m) — (0, 0) (a) m) = (1, 0) (b), (I, m) = (2, 0) 
(c) and (l,m) = (2, 1) (d). The full line is the result of the simulation at T = 0.477. The long dashed line is the result of 
the fixed point iteration for upper cut off l co = 4 at T = 0.600. The short dashed curve is the result of the quasi stationary 
criterion for upper cut off l co = 4 at T = 0.630 
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FIG. 5. Normalized off diagonal elements fo2(l) of the non ergodicity parameter matrix. As normalization Sq Q S 22 was 
chosen. The full line is the result of the simulation at T = 0.477. The long dashed line is the result of the fixed point iteration 
for upper cut off l co = 2 at T = 0.7521. The short dashed curve is the result of the quasi stationary criterion for upper cut off 
I co = 4 at T = 0.630. 
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FIG. 6. Normalized diagonal elements of the matrix of critical amplitudes. As normalization the corresponding 

diagonal elements of the static structure factor are used. The result of the simulation at T = 0.477 is obtained by fitting a von 
Schweidler law plus corrections of order {t/r a ) to the simulation data. The fits to the simulations (full lines) are multiplied 
with a factor 200 (see text for details) to obtain best agreement with the theoretical result for hg (q). Shown are the theoretical 
results for T = 0.7521, upper cut off l co — 2. ^(a)-(c) show the results for (l,m) — (0,0), (1,0), (2,0) respectively. The dashed 
curve is always the result of the theory, ^(d) shows again the same theoretical result as in Fig ^|(c), but the free scale factor 
the critical amplitude of the simulation is chosen to optimize the agreement with theory (see text for details). 
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FIG. 7. Normalized critical amplitude /122(g)- The critical amplitude of the simulation (full line) is multiplied with a free 
scale factor to optimize agreement with theory (see text for details). Shown is the theoretical result for T = 0.7521 (dashed 
line) . 
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FIG. 8. Normalized critical amplitude h% 2 {q)- The critical amplitude of the simulation at T = 0.477 (full line) is multiplied 
with a free scale factor to optimize agreement with theory (see text for details). Shown is the theoretical result for T = 0.7521 
(dashed line). Since the raw data of the simulations had to be smoothed, the data points are omitted. 
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